Quantum noncomutativity in quantum cosmology 
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In the present work, we study the noncommutative version of a quantum cosmology model. The 
model has a Friedmann-Robertson-Walker geometry, the matter content is a radiative perfect fluid 
and the spatial sections have positive constant curvatures. We work in the Schutz's variational 
formalism. We quantize the model and obtain the appropriate Wheeler-DeWitt equation. In this 
model the states are bounded. Therefore, we compute the discrete energy spectrum and the corre- 
sponding eigenfunctions. The energies depend on a noncommutative parameter. We observe that, 
due to the boundary conditions, the noncommutativity forces the universe to start expanding from 
an initial scale factor greater than zero. We also notice that, one can only construct wave-packets if 
the noncommutative parameter is discrete, with a well-defined mathematical expression, in a certain 
region of its domain. 



One important arena where noncommutative (NC) 
ideas may play an important role is cosmology. If su- 
perstrings is the correct theory to unify all the interac- 
tions in nature, it must have played the dominant role 
at very early stages of our Universe. At that time, all 
the canonical variables and corresponding momenta de- 
scribing our Universe should have obeyed a NC algebra. 
Inspired by these ideas some researchers have considered 
such NC models in quantum cosmology [H-Q. It is also 
possible that some residual NC contribution may have 
survived in later stages of our Universe. Based on these 
ideas some researchers have proposed some NC models 
in classical cosmology in order to explain some intriguing 
results observed by WMAP. Such as a running spectral 
index of the scalar fluctuations and an anomalously low 
quadrupole of CMB angular power spectrum An- 
other relevant application of the NC ideas in classical 
cosmology is the attempt to explain the present acceler- 
ated expansion of our Universe una. 

In the present work, we study the noncommutative ver- 
sion of a quantum cosmology model. The model has a 
Friedmann-Robertson-Walker (FRW) geometry, the mat- 
ter content is a radiative perfect fluid and the spatial 
sections have positive constant curvatures. We work in 
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the Schutz's variational formalism (l3| . The noncommu- 
tativity that we are about to propose is not the typi- 
cal noncommutativity between usual spatial coordinates. 
We are describing a FRW model using the Hamiltonian 
formalism, therefore the present model phase space is 
given by the canonical variables and conjugated mo- 
menta: {a, P a ,r, P T }. Then, the noncommutativity, at 
the quantum level, we are about to propose will be be- 
tween these phase space variables. Since these variables 
are functions of the time coordinate t, this procedure 
is a generalization of the typical noncommutativity be- 
tween usual spatial coordinates. The noncommutativity 
between those types of phase space variables have al- 
ready been proposed in the literature. At the quantum 
level in Refs. [JQ, [13] and at the classical level in Refs. 

We quantize the model and obtain the ap- 
propriate Wheeler-DeWitt equation. In this model the 
states are bounded therefore we compute the discrete 
energy spectrum and the corresponding eigenfunctions. 
The energies depend on a noncommutative parameter. 
We observe that, due to the boundary conditions, the 
noncommutativity forces the universe to start expand- 
ing from an initial scale factor greater than zero. We 
also notice that, one can only construct wave-packets if 
the noncommutative parameter is discrete, with a well- 
defined mathematical expression, in a certain region of 
its domain. 



The FRW cosmological models are characterized by 
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the scale factor a(t) and have the following line element, 



ds 2 = -N 2 {t)dt 2 + a 2 {t) 



dr 2 
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where dil 2 is the line element of the two-dimensional 
sphere with unitary radius, N(t) is the lapse function 
and k gives the type of constant curvature of the spatial 
sections. Here, we are considering the case with positive 
curvature k = 1 and we are using the natural unit sys- 
tem, where % = c = G = 1. The matter content of the 
model is represented by a perfect fluid with four- velocity 
= Sq in the comoving coordinate system used. The 
total energy-momentum tensor is given by, 
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where p and p are the energy density and pressure of the 
fluid, respectively. Here, we assume that p = p/3, which 
is the equation of state for radiation. This choice may be 
considered as a first approximation to treat the matter 
content of the early Universe and it was made as a matter 
of simplicity. It is clear that a more complete treatment 
should describe the radiation, present in the primordial 
Universe, in terms of the electromagnetic field. 

Einstein's equations for the metric (fTJ) and the energy 
momentum tensor ^ are equivalent to the Hamilton 
equations generated by the super-hamiltonian constraint, 
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where p a and pr are the momenta canonically conjugated 
to a and T the latter being the canonical variable asso- 
ciated to the fluid [lj]. The total Hamiltonian is given 
by NT-L and we shall work in the conformal gauge, where 
N = a. The commutative version of the present model 
was first treated in Ref. (l5j . 

We wish to quantize the model following the Dirac for- 
malism for quantizing constrained systems (l6j |. First 
we introduce a wave-function which is a function of the 
canonical variables a and T, 



* = *(a,f) 



(4) 



Then, we impose the appropriate commutators between 
the operators a and T and their conjugate momenta P a 
and Pt- Working in the Schrodinger picture, the opera- 
tors a and T are simply multiplication operators, while 
their conjugate momenta are represented by the differen- 
tial operators, 
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Finally, we demand that the operator corresponding to 
NT-L annihilate the wave-function which leads to the 
Wheeler-DeWitt equation, 
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where the new variable r = —T has been introduced. 

The operator NT-L is self-adjoint [l5| with respect to 
the internal product, 



da *(o,r)* $(a, r) , 



(7) 



if the wave functions are restricted to the set of those 
satisfying either ^/(0,r) = or ^'(0, r) = 0, where the 
prime / means the partial derivative with respect to a. 
Here, we consider wave functions satisfying the former 
type of boundary condition and we also demand that 
they vanish when a goes to oo. 

In order to introduce the noncommutativity in the 
present model, we shall follow the prescription used in 
Refs. [1- 3]. In the present model, the noncommutativity 
will be between the two operators a and f, 

[a,f]=ie, (8) 

where a and f are the noncommutative version of the 
operators. This noncommutativity between those opera- 
tors can be taken to functions that depend on the non- 
commutative version of those operators with the aid of 
the Moyal product. Consider two functions of a and 
t, let's say, / and g. Then, the Moyal product be- 
tween those two function is given by: /(a, f) * g(a, f) — 



/(a, r)exp (i0/2)(d&d? - d?d & ) g(a, f) 



Using the Moyal product, we may adopt the following 
Wheeler-DeWitt equation for the noncommutative ver- 
sion of the present model, 
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* vjr(a, t) + 3a 2 * *(a, r) = 0. 



(9) 



It is possible to rewrite the Wheeler-DeWitt equation 
in terms of the commutative version of the opera- 
tors a and f and the ordinary product of functions. In 
order to do that, we must initially introduce the follow- 
ing transformation between the noncommutative and the 
commutative operators, 



a = a - -p f , 



r = r--p & , 



(10) 



and the momenta remain the same. Then, using the 
properties of the Moyal product it is possible to write 
the potential term in Eq. ([9]) in the following way, 



3d 2 *y(d,T) = 3 la- -p? J *(a,f). (11) 

Finally, we may write the commutative version of the 
Wheeler-DeWitt equation ©, to first order in the com- 
mutative parameter 9, in the Schrodinger picture as, 
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For a vanishing 9 this equation reduces to the Schrodinger 
equation of an one dimensional harmonic oscillator re- 
stricted to the positive domain of the variable [HI]. 

In order to solve this equation, we start imposing that 
the wave function ^{a, r) has the following form, 



*(o,r) = A(a)e 



-iEr 



(13) 



Introducing this ansatz in Eq. (|12[) . we obtain the ei- 
genvalue equation, 



d2A ,n 2 „ 
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(14) 



where E is the eigenvalue and it is associated with the 
fluid energy. It is possible to rewrite Eq. (fT4|) such that it 
becomes similar to a one dimensional, quantum mechan- 
ical, harmonic oscillator eigenvalue equation. In order to 
do that one has to perform the following transformations, 
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Introducing these transformations in Eq. ([14]) . we obtain 
the new eigenvalue equation, 



d 2 A 
dx 2 



+ (\-x 2 )A = 0. 



(16) 



This equation is the one dimensional, quantum mechan- 
ical, harmonic oscillator eigenvalue equation and has so- 
lutions for the following discrete values of A, 



A = 2t 



1- 
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where n = 0,1,2,3,.... As a consequency of the second 
transformation of Eq. (|15[) combined with the result from 
Eq. (|17l) . we obtain that the fluid energy is discrete. It 
has the following values, 



E(n,6) = 



36» 2 



-1+ \/l + -(2n+l)6> 2 



(18) 



It is important to notice that E(n, 9) is always positive 
for any value of n and 9, including negative values of 
6. Another important property of E(n, 9) is that, when 
9 -> 0, in Eq. ([18]), we have that E(n,9) -> n + 1/2. 
Which is the correct expression for the one dimensional, 
quantum mechanical, harmonic oscillator energies. For a 
fixed value of 9, E(n, 9) increases when n increases. On 
the other hand, for a fixed value of n, E(n, 9) decreases 
when 9 increases. 

The eigenfunctions to the eigenvalue equation (|16[) are 
given by, 



A n {x) = C n H n {x)e- x2 ' 2 , 



(19) 



where H(x) are the Hermite polynomials of degree n 
and C„ are constans. Since we want to consider solu- 
tions that vanish at the origin, we shall take only the 
odd degree Hermite polynomials: n = 1,3,5,.... Besides 
that, those solutions will be normalized, in the variable 
domain (0,oo), only if the constants C n are equal to: 
C n = 2("- 1 )/ 2 (n!)^ i / 2 (n)- 1 /4. 

We may, now, write the eigenfunctions that are solu- 
tions to the initial eigenvalue equation (Tl4l , with the aid 
of the first transformation of Eq. (|15l) and Eq. (fT5|) . 

A n (a) = C n H n (V6a + 39E/V6)e-^ a+3eE /^ 2 / 2 . 

(20) 

For a fixed n, the wavefunction ^(a, r) Eq. (fT3")) . ob- 
tained with the corresponding eigenfunction A n {a) Eq. 
(f2TJ)) . can only satisfy the boundary condition ^(0, t) = 0, 
when the Hermite polynomial ff„(v6a + 39E/y/&) van- 
ishes. This happens, only, when its argument vanishes, 
which means that, a = —9E/2. Since a and E are posi- 
tive, we learn from that result that 9 has to be negative. 
More importantly, we learn also that in the present non- 
commutative model the scalar factor cannot ever van- 
ishes, a can only vanishes when 9 vanishes, since E is 
always positive. Therefore, the noncommutativity seems 
to prevent the Universe to start from a zero scalar factor 
state. The greater the absolute value of 9, the greater 
the initial value of a. 

The most general expression of ^(a, r) Eq. (|13[) . which 
is a solution to Eq. ([T2"j) , is a linear combination of the 
eigenfunctions A n (a), Eq. (f20|) . combined with the ex- 
ponential factor present in Eq. f)13[) . 
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where the B n 's are constants and E(n, 9) is given in Eq. 
(HHJ). Again, the above wavefunction has to satisfy the 
boundary condition W(0, t) = 0. Observing Eq. (|2~lj) , we 
notice that it will only happen if several Hermite polyno- 
mials, of different degrees n, vanish when their arguments 
vanish. For fixed 9, we notice that for different eigenfunc- 
tions An(a), their arguments vanish for different values of 
a, a n = —9E(n, 9)/2. Therefore, in principle, the ^(a, r) 
Eq. (f2"Tj) cannot satisfy the desired boundary condition. 
In order to obtain a wavefunction that satisfies the above 
boundary condition, we shall make the assumption that, 
in a certain region of the 9 domain, this parameter is 
quantized and is given by: 9 n = a/E(n,6), where a is 
a negative constant. Now, all Hermite polynomials, of 
different degree n, vanish at the same value: a = —a/2. 
Then, ^(a,r) Eq. ([2Tj) satisfies the desired boundary 
condition. 

We may, now, find the explicit expression of 9 n by 
combining the above assumption with the expression of 
E(n,9) Eq. flBJ, 



2a 



(2n + 1) 
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(22) 
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Since, 9 n and a are negative, the denominator of the 
above expression must be positive. From this condi- 
tion , we learn t hat 9 has to be in the following range: 
-^2{2n+ l)/3 < a < 0. As we have seen, above, in 
order for ^(a,r) Eq. (f2"Tj) to satisfy the desired bound- 
ary condition, a has to have a single well defined value. 
We must choose this value, such that, the denominator of 
Eq. (|2"2"j) is always positive, for any value of n. Supposing 
that the eigenfunction of degree one (n = 1) will always 
be present in the linear combination Eq. ([2"T]) . we must 
take the value of a to be in the range: -y/2 < a < 0. 
From Eq. (|22p . we learn that for a fixed value of a, 9 n 
increases when n increases. In the same way, for a fixed 
value of n, 9 n increases when a increases. 

One may understand why the wavefunction must van- 
ishes at a = —a/2 instead of a = by studying the 
potential energy, from the eigenvalue equation (IT41 . If 
we introduce the assumption that 9 n E(n, 9) — a, in Eq. 
(|14|) , we may interpret the following term as the potential 
energy V(a): V(a) = 36a 2 + 36aa. Now, studying this 
potential we learn that it has a minimum not in a = 0, 
but in a = —a/2. The minimum value of V{a) is equal to: 
V(a = —a/2) = —9a 2 . Then, the transformations Eq. 
(|T5|) reduced the present problem with potential V(a) to 
an equivalent harmonic oscillator in the variable x. The 
origin of that harmonic oscillator is located exactly in the 
minimum value of V(a). Therefore, the description of the 
present problem in terms of x must start from x — or, 
in terms of a, must start from a = —a/2. 

Now, we would like to verify that the wavefunction and 
quantities computed with it are well defined, with respect 
to the boundary conditions ^(a = —a/2,r) = ^(a — > 
oo,t) = and the assumption 9 n E(n,9) = a. In order 
to do that, we shall compute, for particular choices, the 
scale factor expected value, (a). First of all, we must in- 
troduce 9 n Eq. (f22|) in E(n, 9) Eq. ([T8]) and obtain an ex- 
pression for the energy as a function of n and a, E(n, a). 
Then, we must define the wavepacket by choosing a fi- 
nite number of eigenfunctions contributing in the linear 
combination Eq. (p?Tj) . Let us choose the first twenty odd 
eingenfunctions, in other words, n = 1, 3, 5, 39, 

39 

*(o,r)= B n A n (a)e- lE ^ T , (23) 

n=l,3,5,... 

Next, we compute the eigenvalues, E(n, a) for those 20 
eigenfunctions, with the aid of Eq. (ITS)) . In order to do 
that we must choose the values of a and the B n . In the 
present situation, we shall choose a = —0.7064965809, 
which is in the range —\[2 < a < 0, B\ = B5 = Bg = 
... = B 37 = 1 and B 3 = B 7 = B n = ... = B w = -1. 



Finally, we must compute the scale factor expected value, 
using the following expression, 
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t 

FIG. 1. (a) (r) for the time interval < r < 50. 

Here, we had to modify the internal product Eq. (|7|), 
due to the fact that the wavefunctions no longer vanish 
at a = 0. Figure 1 shows (a) (r) for the time interval 
< r < 50. It is clear from that figure that (a) (r) 
oscillates between a maximum and a minimum value and 
never vanishes. Therefore, the boundary conditions and 
the assumption we considered give rise to well defined 
wavefunctions and quantities computed with it. 

The fact that (a) (r) never vanishes is expected because 
the integration domain in Eq. (|24l) extends over positive 
values of a and does not include the zero. In the same 
way, if we compute the Ricci and Kretschmann scalars 
expect values [l7|, we would find that they never diverge. 
They will be always finite. In this way, the noncommu- 
tativity acts in a way that prevents, the present model, 
to have an initial singularity at the quantum level. 
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